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Abstract 

Neutrinoless double-beta decay experiments can potentially determine the 
Majorana or Dirac nature of the neutrino, and aid in understanding the neu- 
trino absolute mass scale and hierarchy. Future 76 Ge-based searches target a 
half-life sensitivity of >10 27 y to explore the inverted neutrino mass hierarchy 
Reaching this sensitivity will require a background rate of <1 count tonne -1 
y -1 in a 4-keV-wide spectral region of interest surrounding the Q value of the 
decay. We investigate the overburden required to reach this background goal 
in a tonne-scale experiment with a compact (copper and lead) shield based 
on Monte Carlo calculations of cosmic-ray background rates. We find that, in 
light of the presently large uncertainties in these types of calculations, a site 
with an underground depth >5200 meters water equivalent is required for a 
tonne-scale experiment with a compact shield similar to the planned 40-kg 
Majorana Demonstrator. The required overburden is highly dependent 
on the chosen shielding configuration and could be relaxed significantly if, for 
example, a liquid cryogen and water shield, or an active neutron shield were 
employed. Operation of the Majorana Demonstrator and GERDA de- 
tectors will serve to reduce the uncertainties on cosmic-ray background rates 
and will impact the choice of shielding style and location for a future tonne- 
scale experiment. 

Keywords: neutrinoless double-beta decay, Germanium, cosmic rays 
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1. Introduction 



Neutrinoless double-beta (Qvf3f3) decay may allow the determination of 
the Majorana or Dirac nature of the neutrino, and aid in understanding the 
neutrino absolute mass scale and hierarchy [H El HH IU El E] • Guiding princi- 
ples for extending the Standard Model do not preclude each neutrino mass 
eigenstate from being identical to its antiparticle, or "Majorana" in nature. 
The discovery of Majorana neutrinos would have profound theoretical impli- 
cations for the extension of the Standard Model while yielding insights into 
the origin of mass and supporting the hypothesis that leptogenesis mecha- 
nisms gave rise to the observed baryon asymmetry [7J [8] . 

Neutrino oscillation experiments [HI [101 El] indicate that at least one 
neutrino mass eigenstate has a mass of ~45 meV or more. As a result, 
in the inverted hierarchy mass spectrum with m 3 = meV, the effective 
Majorana mass of the electron neutrino, (nipp), is between ~15 and ~50 meV 
depending on the values of the neutrino mixing parameters and Majorana 
phases. Assuming light Majorana neutrino exchange mediates the decay, the 
limit on the Oz//3/3-decay half-life, T^ 2 , is related to the limit on the effective 
neutrino mass as: 

G» v \M« v \ 2 {m m ) 2 (1) 



1 l/2 



where G 0l/ is an exactly calculable phase-space factor, and Ai 0u is a nuclear 
matrix element. The nuclear matrix element is difficult to calculate, and 
its value is highly model-dependent, but significant progress in reducing the 
theoretical uncertainty has been made recently [12]. In the case that no 
signal is observed, the limit that can be placed on (nipp) can be expressed 
as: 

U.L.(B) 1 1/2 



(2) 



_ln(2)N a TeG°»\M 0u \ 2 

where U.L.(B) is the upper limit on the Ozz/3/3-decay signal given B observed 
background counts, N a is the number of 0v '(3 ,3-decay candidate atoms, T is 
the livetime of the experiment, and e is the signal-detection efficiency. As 
illustrated in Figure [TJ exploring the inverted mass hierarchy by searching for 
Ov(3(3 decay in the isotope 76 Ge requires a sensitivity of Ty 2 > 10 27 y. This is 
a challenging goal requiring tonne-years of exposure and a background rate 
of <1 count tonne -1 year -1 (ct t -1 y -1 ) in the 4 keV-wide spectral region of 
interest (ROI) surrounding the Oz//3/3-decay Q value (Qpp) of 2039 keV, or a 
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rate more than a factor of 100 lower than was achieved in previous-generation 
experiments [TBI H4"] . 

The Majorana PJ Q51 E] and GERDA [H1II9] collaborations are field- 
ing experiments to demonstrate the feasibility of achieving the background 
rates required for the success of a tonne-scale, 76 Ge-based Oi//3/3-decay exper- 
iment. Lessons learned from both experiments will lead to a single, scalable 
design utilizing the best aspects of both technologies. 

The GERDA collaboration plans to operate (in Phase II of their experi- 
ment) a 40-kg array of enriched Ge ( enr Ge) crystals immersed inside a cylin- 
drical liquid argon (LAr) cryostat that is 4 m in diameter and ~5.6 m tall 
|20j . Surrounding the cryostat is a water tank 10 m in diameter and 8.9 m 
tall serving as a neutron moderator and gamma-ray shield. The water tank 
is also a Cherenkov detector equipped with 60 photomultiplier tubes to veto 
events coincident with through-going muons [21J. The collaboration is also 
investigating the possibility of instrumenting the LAr volume to serve as an 
inner active scintillation veto [22J. The background in GERDA is expected 
to be dominated by cosmogenic and primordial radioisotope impurities in the 
Ge crystals and small parts near the detectors |23j . 

The Majorana Demonstrator [T^J HH [17] will be an array of p-type, 
point-contact (PPC) high-purity Ge diodes [HJ [25] housed in two vacuum 
cryostats built from ultra-pure electroformed Cu (EFCu). Each cryostat will 
contain 35 diodes, each with mass of about 0.6 kg. Up to 30 kg of the de- 
tectors will be built from material enriched to at least 86% 76 Ge. A graded, 
passive shield will be constructed from electroformed and commercial high- 
purity Cu, high-purity Pb, polyethylene neutron moderator (some of which 
is borated), and an active plastic scintillator muon veto. The experiment 
will be located in a clean room at the 4850-foot level (4300 meters water 
equivalent, m.w.e.) of the Sanford Underground Laboratory (SUL) at the 
Homestake mine in Lead, South Dakota. The largest sources of background 
in the Majorana Demonstrator are expected to be the decay of cosmo- 
genic isotopes in the Ge, primordial isotope decay originating in the detector 
mounts, cryostat and inner Cu layer of the shield, and inelastic scattering of 
muon-induced neutrons on the Ge, Cu and Pb. 

Both experiments aim to achieve a background rate of 4 ct t _1 y _1 ROD 1 
and either confirm a controversial claim of O^/3/3-decay detection in 76 Ge [26] 
or refute the claim and set a limit on the O^/3/3-decay half-life. In particular, 
the Majorana Demonstrator will set a limit of T®Y 2 > 4 x 10 25 years 
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(90% C.L.) after one year of running. In the case of a tonne-scale experiment 
based on the Majorana Demonstrator concept, additional reduction in 
background is achieved via more powerful analysis cuts associated with larger 
detector arrays, siting the experiment deeper underground, and if necessary, 
using a thicker layer of EFCu in the shield and fabricating the enr Ge detectors 
in an underground laboratory to reduce the amount of cosmogenic activation. 

Here we will focus on cosmic-ray backgrounds and discuss the depth re- 
quirements for a future tonne-scale experiment. In Sec. [2j we describe the 
anticipated sources of background and strategies employed for their mitiga- 
tion. We then quantify the predicted rates of cosmic-ray backgrounds for 
a tonne-scale experiment in the context of the entire background budget, 
drawing upon background studies for a generic 0^/3,3-decay experiment with 
a "compact" shield similar to the Majorana Demonstrator, in Sec. [3} 
In Sec. [4] we derive the depth requirements for a tonne-scale experiment and 
discuss the interplay between the required depth of the laboratory and the 
style of shielding. 

2. Backgrounds &z Mitigation Strategies 

Achieving a low background is critical to the search for 0v(3(3 decay. In 
76 Ge, 0v(3(3 decay would produce a mono-energetic peak at 2039 keV. With 
a typical Ge detector resolution of 0.2% FWHM at 2 MeV, less than a frac- 
tion of ~2 x 10~ 15 of the continuous two-neutrino double-beta {2vf3f3) decay 
spectrum contaminates the ROI, making this class of background completely 
negligible [2]. Any particle with E > Qpp originating outside of the Ge 
detector volume depositing a portion of its energy presents a potential back- 
ground. Even though the excellent inherent energy resolution of Ge detectors 
serves to reduce the count rate in the ROI from many sources of background, 
the detectors themselves and the materials used to fabricate the detector 
mounts, electronic readout, and shielding must be of extreme radiopurity. 
The experiment must employ active and passive shielding to reduce external 
backgrounds originating from natural radioactivity in the cavern walls, and 
must be sited underground to reduce cosmic-ray backgrounds. 

2.1. Depth- dependent backgrounds 

Cosmic-ray muons produce electromagnetic (EM) and hadronic particle 
cascades as they pass through rock, shielding materials, and the Ge diodes. 
Highly-penetrating neutrons with energies extending up to several GeV can 
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be generated in these cascades. In the Majorana Demonstrator, these 
neutrons can produce backgrounds via elastic scattering on Ge, Ge(n, n), 
and inelastic scattering, (n, n/7), on Ge, Pb and Cu. In addition, about 
a dozen unstable isotopes of Ge, Ga and Zn with half-lives greater than a 
few seconds and with Q > Qpp can be produced in the Ge diodes in situ 
by spallation processes. Rarely, muons stop inside the detector and form 
muonic atoms in which the muon cascades down into the IS state and either 
decays via — > e'VeV^ or is captured by the atomic nucleus via the weak 
interaction process pT + p — > n + |27J . The reaction typically results in 
nuclear excitations of ~100 MeV followed by the emission of gamma rays 
and a number of neutrons with energies typically less than 10 MeV but as 
high as a few tens of MeV as the daughter nucleus de-excites. 

Neutron inelastic scattering interactions producing gamma rays coinci- 
dent with muon passage through an anti-cosmic veto are relatively easy to 
eliminate. Of greater concern are un-vetoed neutrons resulting from muon 
interactions in the walls and other components of the experimental cham- 
ber. While the /z-induced neutron flux is typically at least three orders of 
magnitude lower than the neutron flux from primordial radioactivity at un- 
derground depths greater than ~3000 m.w.e. [28], these neutrons are more 
penetrating and thus more difficult to shield against. Muon- induced neutron 
backgrounds can be suppressed by choosing a site deeper underground, or by 
minimizing the use of high-Z materials in the shield. 

As depth increases, the muon energy spectrum becomes harder and the 
integral flux decreases. As a rule of thumb, the muon flux decreases by 
about an order of magnitude for every 1500 m.w.e. increase in depth [28J. 
The integral muon flux at depths greater than ~1 kilometer water equivalent 
(km.w.e.), 0^, can be parameterized as 

( p iM (h) = a e~ h ^ 285 + a 1 e- h /° m8 (3) 

where h is the flat-overburden-equivalent depth in km.w.e., = 68 x 
10 -6 cm -2 s _1 and a\ = 2.1 x 10 _6 cm -2 s^ 1 [28]. Since the muon energy 
spectrum becomes harder with increasing depth, the /i-induced neutron flux 
decreases somewhat more slowly with increasing depth and can be parame- 
terized as 

Mh) = Po(^)e- h/Pl (4) 
where P = (4.0 ± 1.1) x 10~ 7 cm" 2 s" 1 and P 1 = 0.86 ± 0.05 km.w.e. [28]. 
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Furthermore, since the cross section for the in-situ production of cosmogenic 
isotopes increases as a positive power of the muon energy, the rate of cos- 
mogenic isotope production also generally does not decrease as fast as the \i 
flux with depth. Several studies of //-induced isotopes in materials ranging 
from liquid scintillator to Pb have found that the production cross section, 
a, goes as 

a(E,) oc El (5) 

where is the muon energy and a is a constant found to range between 
0.50 and 0.93 with a mean value (a) = 0.73 ± 0.10 [29]. Using this relation, 
one can show [2E] that the rate of /i-induced cosmogenic isotope production, 
Riso, scales with depth as 

Risojh = 0) = / 4GeV y ^(fr = 0) 
Riso(h) V (E,) J Mh) 

where (E^) is the mean muon energy at depth h, which can be written as 

The parameters e^, b and 7^ depend on the muon energy spectrum in the 
atmosphere, muon energy loss processes in rock, and local rock density and 
composition. Values provided by Lipari & Stanev [30] (b = 0.383 / km.w.e., 
7 M = 3.7 and e M = 618 GeV) and Battistoni et al. [31J (b = 0.4 / km.w.e., 
7 M = 3.77 and e M = 693 GeV) both agree with the measured value of the 
mean muon energy at the Gran Sasso National Laboratory (LNGS) [32J. 

2.2. External & intrinsic backgrounds 

Significant backgrounds arise from naturally-occurring primordial radioac- 
tivity present in materials, including in the walls of the underground cavern. 
One of the isotopes in the 232 Th decay chain, 208 T1, is of particular concern, 
since it decays with the emission of a 2615-keV gamma ray, well above Qpp 
for 76 Ge. An isotope in the 238 U decay chain, 214 Bi, also decays with the 
emission of gamma rays with energies above Qpp. Radon and the plate-out 
of its daughters inside the shield and cryostat present a variety of potential 
a, (3 and 7 backgrounds. A massive bulk shield composed of radiologically 
clean material and a Rn exclusion box surrounding the detector systems 
has been repeatedly demonstrated to greatly reduce primordial gamma-ray 
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backgrounds originating externally to the shield. 

Neutrons resulting from the decay of primordial isotopes in the cavern 
walls, such as from the spontaneous fission of 238 U or from natural alpha 
emitters via (a, n) reactions, are also of concern. These neutrons can pro- 
duce gamma rays with energies above Qpp through inelastic scattering and 
capture processes in the shielding materials or in the detectors themselves. 
With a maximum energy of ~10 MeV, these neutrons can be substantially 
suppressed by incorporating a hydrogenous neutron absorber into the bulk 
shield. 

Great care must also be taken to minimize the radioactive contaminants 
present in the components used inside the shield. Copper has become a 
material of choice for the construction of ultra-low-background detectors be- 
cause of its lack of naturally-occurring radioisotopes, its excellent physical 
properties, and the ability to purify it to a high degree via electrodeposition 
[33| [Ml [35] . One must ensure that it does not have significant quantities 
of cosmogenic radioisotopes, the most troublesome of which is 60 Co, gener- 
ated by (n, a) reactions on 63 Cu. The strategy employed by Majorana to 
mitigate this background is to electroform Cu parts underground. 

Germanium detectors have the advantage that problematic radioisotopes 
are effectively removed during the detector production process. The cos- 
mogenic isotope 68 Ge is removed during the enrichment of the Ge material. 
Impurities of U and Th, as well as cosmogenic 60 Co, are effectively removed 
during the zone refinement and crystal pulling stages. The relatively long- 
lived cosmogenic isotopes 60 Co (Ti/2 = 5.3 y) and 68 Ge (Ty 2 = 271 d) are 
produced via cosmic-ray spallation in Ge. With a Q value of 2.8 MeV, 60 Co 
j3 decay contributes to the background continuum. Although 68 Ge decays by 
electron capture with Q < Qpp, its daughter isotope, 68 Ga, has a half-life 
of 68 minutes and decays predominantly by /3 + emission with a Q value of 
2.9 MeV. The K- and L-shell X rays from the 68 Ge decay can potentially 
be used to tag the subsequent 68 Ga decay and veto any events occurring 
within the time span of a few 68 Ga half-lives. In addition, 68 Ge has a half-life 
much shorter than the expected ~ 10- year time scale of a tonne-scale exper- 
iment. More direct mitigation of cosmogenic backgrounds must be achieved 
by carefully limiting the aboveground exposure of the enriched material and 
completed detectors to cosmic rays, or (more ideally) by conducting as much 
of the detector fabrication process as practicable in an underground labora- 
tory [36J. 
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3. Simulations & Background Studies 



In this section we synthesize the results from a number of simulations and 
background studies to estimate the total background for a tonne-scale Ovf3{3- 
decay experiment with a compact shield. The rate of cosmic-ray backgrounds 
as a function of depth for 76 Ge-based Ov (3 /3-decay experiments was studied 
by Mei & Hime [28] with an extensive FLUKA- and GEANT3-based custom 



Monte Carlo (MC) calculation. In Sec. 3.1 we summarize the salient results 



and use the scaling relations in Sec. 2.1 to derive cosmic-ray background rates 



for an experiment situated at the 4850-foot level of SUL. Additional results 
attained using the same framework but performed after the publication of 
Mei & Hime [28] are used to describe how modifications in the thickness of 
neutron shielding and muon veto efficiency change the results. In Sec. 32 we 
describe our simulations of a 60-kg detector module with a compact shield 
similar to the Majorana Demonstrator conducted within the MaGe [3T] 
simulation framework. Based on the simulation, we use a parametric model 
to link material radiopurity goals for various shield and detector components 
to background counts in O^/3/3-decay ROI, taking into account analysis cut 
efficiencies, to demonstrate that a rate of ~1 ct t _1 y _1 ROI -1 can be achieved 
in a tonne-scale experiment. 

3.1. Cosmic-ray muon-induced backgrounds 

Simulated cosmic-ray background rates, including direct contributions 
from untagged muons, muon capture, untagged //-induced fast neutron elas- 
tic and inelastic scattering, and cosmogenic radioactive isotope production in 
situ, have been reported by Mei & Hime [2S] . Their work employed FLUKA 
to simulate the interactions of cosmic ray muons in the rock surrounding the 
experimental apparatus. The muon energies and angular distributions were 
sampled from empirical parameterizations. The generated spallation neu- 
trons were propagated to the cavern walls, at which point they were booked 
in histograms to be used in event generators for a GEANT3-based simulation 
of the detector response, using Hauser-Feshback theory to calculate neutron 
inelastic scattering cross sections. The simulated experimental set-up in- 
cluded a 60-kg array of enriched Ge crystals inside a compact shield with 
a 10-cm-thick inner Cu layer followed by 40 cm of Pb, 10 cm of polyethy- 
lene, and an outer active muon veto with 90% efficiency. To enable an easy 
comparison to other experiments, the simulated set-up was situated at the 
depth of LNGS (3100 m.w.e., flat-overburden-equivalent). Mei & Hime [28] 
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found that the dominant //-induced contribution to the background is from 
fast neutron elastic and inelastic scattering on Ge. The background rates of 
the simulated processes are summarized in Table [T] 

Before discussing how the results scale with changes in the shield con- 
figuration, it is necessary at the outset to mention that MC estimates of 
/i-induced background rates at great depths underground are fraught with 
significant uncertainty. An uncertainty of 10-15% on the muon flux at depth 
is attributed to the uncertainties in our knowledge of primary muon flux 
and energy spectrum in the atmosphere [28, 38J. Uncertainty in the effec- 
tive depth of the underground laboratory results in 30% uncertainty in the 
muon flux normalization for the case of the 4850-foot level of SUL (4300 
± 200 m.w.e. [28]). Kudryavtsev et al. [38] compared differences in muon 
propagation for two different types of rock composition ( "standard rock" and 
Modane rock) and found a 15% change in the muon flux normalization at a 
depth of 5000 m.w.e., but no significant change in the mean muon energy. 
Adding these uncertainties in quadrature results in a 37% uncertainty on the 
muon flux normalization at 4300 m.w.e. depth. 

Uncertainties on the //-induced neutron fluxes are larger. The /i-induced 
neutron yield in materials with Z < 40 disagrees by 40% between GEANT4 
and FLUKA and by a factor of 2 for high-Z materials such as Pb [39]. While 
neutron yields in FLUKA are higher than in GEANT4, comparison to ex- 
perimental data indicates that both codes under predict neutron yields in 
light materials, and probably do so for Pb as well [39]. We therefore include 
a factor-of-two uncertainty on the /i-induced neutron yield when calculating 
our background rates. Uncertainties due to neutron tracking in the Monte 
Carlo are typically taken to be ~20% based on the level of agreement between 
different codes [40] ■ Mei & Hime [2H] note that neutron backscattering off of 
the walls of the experimental cavern can result in a factor of 2 to 3 change in 
the neutron flux incident on the shielding depending on the exact dimensions 
of the cavern and shield. Mei & Hime [28] included the effects of neutron 
backscattering in their simulations for a cavern with dimensions of 30 x 6.5 
x 4.5 m 3 . We take the uncertainty on the neutron flux due to the effect of 
backscattering to be 150%, and note that for a larger cavern, the contribution 
of backscattering to the total muon-induced background is expected to be 
smaller. Adding the logarithms of the uncertainties in quadrature, we find a 
total uncertainty on the muon-induced neutron flux of a factor of 3.2. 

With these limitations in mind, we scale the Mei & Hime [2H] results to 
reflect the /i-induced background rates in a more realistic representation of 
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the planned Majorana Demonstrator shield configuration. The major 
differences are: a thicker (30-cm) layer of polyethylene neutron moderator, 
a 99%-efficient muon veto, and an overburden corresponding to the depth of 
the 4850-foot level at SUL (4300 ± 200 m.w.e.). The fraction of unvetoed 
direct muon hits is decreased by a factor of 10 on account of the more ef- 
ficient anti-cosmic veto. According to Eqs. [3] and [4] we also gain reduction 
factors of 5.9 and 5.6 for the direct muon flux and //-induced neutron flux, 
respectively, on account of the greater depth. Since the "Others" category 
in Table Q] includes in- situ production of cosmogenic isotopes in addition to 
some processes whose rates decrease more quickly with depth, we conserva- 
tively use Eq. [6] to quantify the depth dependence of backgrounds in that 
category, resulting in a reduction factor of 5.3. We have assumed that the 
effects of changes in the energy spectrum of neutrons penetrating through 
the thicker layer of polyethylene can be neglected in comparison to the large 
uncertainties on the overall /i-induced neutron flux. Simulations performed 
after the publication of Mei & Hime [2S] using the same framework showed 
that increasing the thickness of the polyethylene from 10 cm to 30 cm re- 
duces the //-induced fast neutron flux penetrating the shield by a factor of 
two, and that the more efficient anti-cosmic shield decreases the number of 
untagged neutrons from muons that do not cross the detector by a factor of 
1.8. Even though most /i-induced neutrons are accompanied by the muon 
itself, the detection efficiency is much lower for those that are not, and so 
we also conservatively use this factor to describe the improvement in veto 
efficiency for all /t-induced neutrons. 

Muon-induced neutrons emerging from the cavern wall are typically ac- 
companied by electromagnetic (EM) cascades. The Mei & Hime (2E] simula- 
tion treated the generation of /i-induced neutrons in the rock separately from 
their propagation into and through the detector shielding. Simulations of the 
EDELWEISS experiment, a Ge-based dark matter experiment with a com- 
pact shield design similar to that envisioned here, indicate that including EM 
cascades improves the veto efficiency by almost an order of magnitude |4"T] . 
Given the differences in the backgrounds relevant for dark matter searches 
and 0i>PP decay, we use this value as an upper limit and conservatively choose 
a factor of four background reduction, with a factor-of-two uncertainty, to 
correct for the missing EM cascades in the Mei & Hime (28] simulations. The 
corresponding uncertainty introduced on the background rates is 60%, much 
smaller than the uncertainty due to neutron yield and tracking. By including 
this factor, any depth-dependent change in neutron veto efficiency due to the 
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^-induced neutron energy spectrum becoming harder with increasing depth 
does not significantly change the results. The scaling factors used to modify 
the Mei & Hime [28] results as discussed above are summarized in Table |2j 
Mei & Hime (28] state that the application of analysis cuts can poten- 
tially reduce the total /z-induced neutron background rates by a factor of 
~7.4. However, this factor should decrease as the anti-cosmic veto efficiency 
increases from 90%, as assumed in Mei & Hime |28j, to 99% as assumed in 
this work. The granularity and pulse-shape-analysis (PSA) cuts (see Sec. 



3.2) are most effective for rejecting high-multiplicity events, but as the veto 
efficiency increases, the average multiplicity of events that do not trigger the 
veto decreases. Based on our experience simulating other backgrounds, we 
expect analysis cuts to contribute at least a factor 2.5 reduction (cut effi- 
ciency of 40%) in the rate of /i-induced neutron backgrounds. Assuming this 
corresponds to an upper limit on the cut efficiency, and assuming the factor 
of 7.4 from Mei Sz Hime [28] corresponds to a lower limit (14% cut efficiency), 
we will assume an analysis cut efficiency of 25%, with an uncertainty of 60%. 

Including the effects of analysis cuts and logarithmically adding in quadra- 
ture the uncertainties on the cut efficiencies (60%), muon flux normalization 
(37%), neutron yield and tracking (a factor of 3.2), and the scaling factors 
introduced to adjust the Mei & Hime [28J results (60%), the total /i-induced 
background rate for an experiment with a compact shield configuration situ- 
ated at 4300 m.w.e. is 0.33 ct t _1 y _1 ROR 1 with an uncertainty of a factor 
of ~4. 

3.2. MaGe simulations 

MaGe is an object-oriented simulation package based on ROOT [42J and 
the GEANT4 toolkit [43] 144] optimized for simulations of low-background 
Ge detector arrays [37]. MaGe is developed and maintained jointly by the 
Majorana and GERDA collaborations. Simulations performed to repro- 
duce experimental data taken with a number of different detectors and ra- 
dioactive sources generally reproduce spectral features - such as relative peak 
heights, peak to Compton ratios, and Compton scattering multiplicity - to 
within 5-10% [37] . In cases where GEANT4 has been shown to perform 
inadequately when compared to data, such as in neutron production from 
cosmic-ray muons and neutron inelastic scattering [45] , efforts are underway 
within the Majorana, GERDA, and GEANT collaborations to correct or 
compensate for the inadequacies. 
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Neutrinoless double-beta decay results in the emission of two electrons, 
and because these electrons travel very short distances (~1 mm) in Ge, their 
energy depositions are effectively single-site events. In contrast, most back- 
grounds arise from the interaction of gamma rays in the detector. Gamma 
rays with energies at or above Qpp, with a mean distance of centimeters be- 
tween Compton scatters in Ge, are likely to deposit only a portion of their 
energy. If the energy deposited is within a few keV of Qpp, the interac- 
tion represents a background to the Oz//3/3-decay signal. Such interactions 
are typically multi-site energy depositions. Multi-site events can be identi- 
fied and rejected if the multiple scatters occur in more than one crystal in a 
tightly-packed array. This granularity technique is estimated to contribute a 
factor 2 to 5 suppression of multi-site events, depending on the origin of the 
background. Modern signal digitization techniques also allow for the discrim- 
ination between single-site and multi-site event topologies in a single crystal. 
P-type, point-contact detectors are particularly suitable for the application 
of these techniques. Measurements indicate that PSA techniques applied to 
data from PPC detectors can reject >90% of multi-site background events 
while retaining >90% of single-site events [161 Wf\. 

We have employed the MaGe simulation package to determine the back- 
ground rates from each major background source discussed in Sec. [2] for 
generic 76 Ge-based Oz//3/3-decay experiments. The simulated geometries in- 
cluded 3-, 9-, 21- and 57-detector arrays of 1.05-kg, semi-coaxial Ge crystals, 
and a 35-crystal array of 560-g PPC detectors. For each geometry, the Cu 
mounts, front-end read-out electronics near the detectors, and cables were 
simulated inside a Cu cryostat. The detector module was enclosed inside a 
compact shield of Cu and Pb corresponding to the baseline thicknesses for 
the Majorana Demonstrator described above. Using MC calculations 
of the detector response, measured or targeted levels of impurities can be 
linked to background rates in the 4-keV-wide ROI. For example, in the 232 Th 
and 238 U decay chains, we focus on the two radionuclides that have gamma- 
ray lines above Qpp: 208 T1 and 214 Bi. The estimates include component- and 
background-specific suppression factors, including those from the granularity 
and PSA cuts, and for 68 Ga decays, a time-correlation cut corresponding to 
five 68 Ge half lives. The granularity cut is applied to any event in which 
more than one detector is triggered and records >5 keV of energy deposited. 
Since our PSA software was not fully implemented at the time when the 
simulations were run, the PSA background suppression capabilities were esti- 
mated by applying a heuristic based on the charge-drift properties and signal 
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characteristics of PPC detectors, including the impact of noise and readout 
bandwidth. PSA simulations were subsequently performed and the results 
were found to be consistent with the heuristic-based analysis. The results 
for the 57-crystal array (the size envisioned for a tonne-scale experiment) are 
summarized in Tabled) 

Based on the absence of clear a peaks in the Heidelberg-Moscow data 
|13j . we place upper limits on the 238 U and 232 Th contamination in the Ge 
crystals of <1 x 10 -15 g/g and <4 x 10 -15 g/g, respectively, corresponding to 
a background rate of <0.30 ct t -1 y -1 ROI 1 from all decay processes in the 
two decay chains, assuming secular equilibrium. We estimate a background 
rate of <0.04 ct t -1 y -1 ROI -1 from the decay of cosmogenic 68 Ge and 60 Co 
in the Ge crystals assuming the enr Ge material is limited to a maximum of 
100 d of sea-level-equivalent cosmic-ray exposure. Our estimate includes the 
effects of PSA and granularity cuts, and the veto of any events occuring 
within five 68 Ga half-lives after the observation of a K-shell X ray from the 
decay of 68 Ge. The largest contribution to the total background (0.46 ct t -1 
y -1 ROI -1 ) is expected to originate from the EFCu portion of the shield and 
the cryostats. Additional sources of background present inside the shield 
include the detector mounts, front-end electronics, cables, and other small 
parts, contributing <0.14 ct t -1 y -1 ROI -1 . 

By simulating different background sources under a variety of detector 
shielding configurations, we were able to parametrically extrapolate back- 
ground estimates for different shielding designs. By doubling the thickness of 
EFCu shielding relative to that proposed for the Majorana Demonstrator, 
we estimate that backgrounds in a tonne-scale experiment from an outer, 
commercial Cu shielding layer would be ~0.02 ct t -1 y -1 ROI -1 . In this 
configuration, backgrounds from natural radioactivity in the Pb would be 
negligible (<0.02 ct t -1 y -1 ROI -1 ). The effectiveness of the shield and 
neutron moderator limits the background from external gamma rays and 
neutrons from the cavern walls to ~0.05 ct t -1 y -1 ROI -1 . 

The total depth-independent background rate of <1.1 ct t -1 y -1 ROI -1 
carries an uncertainty of 30%. We note that the sum of depth-independent 
backgrounds leaves no room for depth-dependent backgrounds if the back- 
ground goal of <1 ct t -1 y -1 ROI -1 is to be achieved. 
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4. Depth Requirements 

A tonne-scale, 76 Ge-based Oz//3/3-decay experiment must strive to make 
//-induced backgrounds negligible compared to the ~1 ct t _1 y -1 ROR 1 ex- 
pected from depth- independent backgrounds. At a depth of 4300 m.w.e., the 
rate of /i-induced neutron backgrounds for a compact shield configuration 
similar to the Majorana Demonstrator is expected to be ~0.3 ct t _1 
y -1 ROR 1 after analysis cuts, but with an uncertainty of a factor of 4, domi- 
nated by uncertainties associated with the //-induced neutron yield in high-Z 
materials and the effects of neutron backscattering. A target /i-induced back- 
ground rate of 0.1 ct t _1 y _1 ROR 1 implies a depth requirement of >5200 
m.w.e. assuming the expectation value of the estimated background rate. 
However, given the large uncertainty, an even deeper location would be pre- 
ferred to mitigate the associated risk. Figure [2] illustrates this requirement in 
the context of several existing and possible future underground laboratories. 

Table [T] indicates that the depth-dependent background for an experiment 
with a compact, high-Z shield configuration is dominated by /i-induced neu- 
tron inelastic scattering on Cu and Pb components of the shield. Eliminating 
high-Z materials close to the Ge detectors reduces this background contribu- 
tion and can significantly relax the depth requirement. This is the approach 
taken by GERDA. The rates of /i-induced backgrounds in GERDA at the 
depth of LNGS (3100 m.w.e.) have been calculated using a MC simulation 
by Pandola et al. [40]. The production and subsequent decay of cosmogenic 
isotopes in the Ge crystals and the LAr volume was found to be the limiting 
depth-dependent background with a rate of ~0.8 ct t _1 y _1 ROR 1 , including 
a conservative accounting for the effects of analysis cuts. The decay of 77 Ge 
and 77m Ge produced by neutron capture of thermalized /i-induced neutrons 
produced in the shield volume was found to be the dominant contribution. 
The authors assign statistical and systematic uncertainties of ~20% and 
45%, respectively, to the estimated rate. The lower systematic error com- 
pared to our estimate for a compact shield is possible because for materials 
with Z < 40, the predicted neutron yields are in closer agreement between 
different MC codes than those for Pb, and backscattering neutrons in the 
cavern are not expected to contribute significantly to the background. Pan- 
dola et al. [40] also found that the /t-induced background rates are an order 
of magnitude smaller in the case that liquid nitrogen is used in the shield in 
place of LAr [40] because the /i-induced particle showers do not develop to 
the same extent as in LAr. 
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The /i-induced background rate in GERDA would be reduced by a factor 
~5 to ~0.15 ct t" 1 y" 1 ROR 1 at the depth of SUL according to Eq.g When 
comparing this rate to our estimate of 0.33 ct t _1 y _1 ROR 1 for a compact 
shield, we must keep in mind that the uncertainty is significantly smaller, 
and that the estimated rate for GERDA is an upper limit since the assumed 
analysis cut efficiency for the 77 Ge background is probably conservative [40J . 

The depth requirement for an experiment with a compact shield configu- 
ration is driven in large part by the present uncertainty associated with MC 
predictions of neutron yields in high-Z materials. The depth requirement 
can be relaxed significantly for a low-Z shield configuration since the asso- 
ciated /i-induced background rates for low-Z configurations are lower and 
have smaller uncertainties. In particular, for a GERDA-style experiment 
using LN 2 for shielding, //-induced backgrounds would be negligible at the 
depth of the SUL 4850-foot level. It should be noted, however, that such a 
shield configuration would be more technically complex, would require much 
more space, and could present significant oxygen deficiency and other safety 
hazards. 

5. Conclusions 

The Majorana and GERDA collaborations are developing experiments 
based on different shielding techniques to demonstrate that a tonne-scale, 
76 Ge-based Oi//5/5-decay experiment can achieve the extremely low background 
rates necessary to accomplish the science goals. Achieving the desired scien- 
tific impact depends crucially upon reducing the background rate to <1 ct 
t _1 y _1 ROD 1 . Even with the utmost attention to cleanliness and the ra- 
diopurity of construction materials, a tonne-scale experiment will likely have 
to contend with a background rate of ~1 ct t _1 y _1 ROR 1 from impurities 
in the detectors, their mounts and cables, and the shield. Therefore, a tonne- 
scale experiment should strive to make negligible all /i-induced backgrounds. 

The planned depth of 4300 m.w.e. for the Majorana Demonstrator is 
ideal in that /z-induced backgrounds are expected to be measureable but not 
dominant, thus enabling a proper assessment of backgrounds from various 
different sources. We have shown that /i-induced backgrounds are expected 
to be of significant concern for a tonne-scale experiment with a shield similar 
to that of the Majorana Demonstrator if it were operated at the same 
depth. 
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A tonne-scale Oz//3/3-decay experiment with a compact shield will likely 
need to be located at a depth of >5200 m.w.e. Based on our current knowl- 
edge, an even deeper location is preferred to mitigate the risks arising from 
large uncertainties in the estimated /i-induced background rates stemming 
in large part from poor knowledge of /^-induced neutron yields in high-Z ma- 
terials. We note that the largest single uncertainty in the background rate 
corresponds to neutron backscattering whose effects are highly dependent on 
the cavern and experimental geometry. For a cavern larger than that simu- 
lated in Mei & Hime [28], this uncertainty should contribute to a downward 
change in the estimated muon-induced background rate. The depth require- 
ment could also be significantly relaxed if a low-Z shielding configuration or 
an active neutron veto system were employed. 

The eventual choice of the location for a tonne-scale experiment will 
involve a careful analysis of the backgrounds achieved in the Majorana 
Demonstrator and GERDA together with an evaluation of the complex 
interplays between the cosmic-ray shielding afforded by depth, technical fea- 
sibility, and background suppression capability of alternative shield designs. 
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Figure 1: The sensitivity of a search for 0vj3j3 decay using 86%-enriched 76 Ge as a function 
of background rate and exposure calculated according to equations [l] and [2] We take the 
signal detection efficiency, e, to be 68% to account for the effects of analysis cuts. The 
exposure is defined to as the product of the active mass of the detectors and the live time 
of the experiment. The phase space factor G 0v is taken from Ref. [IS] and the conversion 
from half-life to (mpp) is performed using a QRPA-calculated nuclear matrix element 
[49] . The dark-green band indicates the range of {mpp) in the inverted mass hierarchy 
for the full range of CP-violating phases, while the light-green bands take into account 
uncertainties in the measured neutrino oscillation parameters. The half-life corresponding 
to the recent claim of O^/3/3-decay detection [33] is shown in blue. 
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Figure 2: The //-induced background rate for a tonne-scale O^/3/3-decay experiment with 
a compact shield is plotted according to Eqj3l normalized to a rate of 0.33 ct t _1 y _1 
ROI -1 at 4300 m.w.e. (see Table [l] and Sec. 3.1). The blue band represents a factor of 4 
uncertainty in the background rate. The diagonally-shaded bands show the effective depths 
(and associated uncertainties according to Ref. [28]) for several underground laboratories. 
The target /i-induced background rate corresponds to 0.1 ct t^ 1 y _1 ROI -1 . 
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Reaction 


BG in Ref. [28J 


Scaled BG 


76 Ge(n, n'7) 


40 


0.49 


74 Ge(n, n'7) 


8.0 


0.10 


Cu(n, n'7) 


7.6 


0.094 


208 Pb(n,n' 7 ) 


14 


0.17 


Ge(n, n) 


14 


0.17 


\i hits 


10 


0.17 


Others 


9.6 


0.13 


Total 


100 


1.3 



Table 1: Muon-induced background rates for a generic O^/3/3-decay experiment with a 
compact shield similar to the planned Majorana Demonstrator operating at Gran 
Sasso (3100 m.w.e.) reported in Mei & Hime [33] and the same background rates scaled as 
described in Sec. [XT] to represent a more realistic detector configuration at the 4850-foot 
level (4300 m.w.e.) of Sanford Underground Laboratory. The "Others" category includes 
muon capture, neutron capture on Ge and Cu, and prompt backgrounds due to cosmogenic 
isotopes produced in situ. Rates are given in ct t _1 y _1 ROI -1 and do not include the 
effects of analysis cuts. 



Background Depth n Moderator fi/n Veto Total 



*(n, n'7) 


5.6 


2 


7.2 81 


Ge(n, n) 


5.6 


2 


7.2 81 


fj, hits 


5.9 




10 59 


Others 


5.3 


2 


7.2 76 



Table 2: Summary of the scaling factors (described in Sec. 3.1 ) used to derive the back- 
ground rates shown in the right-hand column of Tabic [TJ The scaling factors listed for 
*(n, n'j) apply to all neutron inelastic scattering processes in Table [TJ The factors of 7.2 
reflect the combined effects of a higher veto efficiency (a factor of 1.8) and accounting for 
the missing EM cascades accompanying muon-induced neutrons emerging from the cavern 
walls in the simulations (a factor of 4). 
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Source 


Radioactive Isotope 

[ct t- 1 y- 1 ROI- 1 ] 


Background 

[ct t- 1 y- 1 ROI" 


onr Ge Crystals 


««Ge 


BU Co 


232 Th / 

238|j 


<0.34 




<0.01 


<0.03 


<0.30 




Detector Mounts 


^riQ m i 

zUo r pi 

0.02 


214 B j 

0.04 




0.06 


Front-end Electronics 


208 rpi 

1 1 
0.01 


214™; 

r>i 
0.02 




0.03 


Cables 


208rpj 

0.02 


214 Bi 

0.03 


60 Co 
<0.01 


<0.05 


Cryostat 


208 T j 

0.12 


214 Bi 

0.03 




0.15 


Inner Cu Shield 


208 T j 

0.24 


214 Bi 

0.07 




0.31 


Outer Cu Shield 


208rpj 

0.01 


214 Bi 

0.01 


60 Co 
<0.01 


<0.03 


Pb Shield 


208 T j 

<0.01 


214 Bi 

<0.01 




<0.02 


Other 


Surface a 


Rock 


Other 


0.11 


0.05 


- 0.05 


- 0.01 


Total Depth-independent 








<1-1 


Total Depth-dependent 








-0.02 


Total Background 








<1-1 



Table 3: Background budget for a tonne-scale detector based on simulations of a 60- 
kg, 57-crystal detector module in a compact shield. The depth-dependent background 
assumes the experiment is located at the 7400-foot level (—6600 m.w.e) of the Sanford 
Underground Laboratory at Homestake. The total depth-independent background rate 
carries an uncertainty of 30%. 



25 



